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We report on a recent investigation of the single spin asymmetry (SSA) in low virtuality 
electroproduction of J/rp in color evaporation model. We show that this can be used as a 
probe for the still unknown gluon Sivers function. 



1 Introduction 

Large single spin asymmetries (SSA) observed when an unpolarized beam of electrons or pro- 
tons is scattered off a transversely polarized target can be explained with the inclusion of k± 
dependence in parton distribution functions (pdf's) and fragmentation functions (fPs) Q]. One 
is led to a generalized factorization formula called TMD factorization [3J [3], which in some 
processes has been proved at leading twist and leading order [4] and has been argued to hold at 
all orders. The inclusion of the effect of transverse momentum of partons in parton distribution 
(pdf) and fragmentation functions leads to a new class of parton distributions that include the 
effects of spin and transverse momentum of the partons. One of these functions is the Sivers 
function which describes the probability of finding an unpolarized parton inside a transversely 
polarized hadron. In this work, we propose charmonium production as a probe to investigate 
the Sivers function and as a first step, estimate SSA in photoproduction (low virtuality elec- 
troproduction) of charmonium in scattering of electrons off transversely polarized protons. At 
leading order (LO), this receives contribution only from a single partonic subprocess jg —> cc . 
Hence, SSA in e + — >• e + J/tp + X, if observed, can be used as a clean probe of gluon Sivers 
function. In addition, charmonium production mechanism can also have implications for this 
SSA and therefore, its study can help probe the production mechanism for charmonium. 

2 Estimate of the Sivers Asymmetry 

There are several models for charmonium production. We use the color evaporation model 
(CEM) as its simplicity makes it suitable for an initial study of SSA in the charmonium pro- 
duction. This model was first proposed by Halzen and Matsuda [5] and Fritsch [6]. In this 
model, a statistical treatment of color is made and the probability of finding a specific quarko- 
nium state is assumed to be independent of the color of heavy quark pair. In later versions 
of this model it has been found that the data are better fitted if a phenomcnological factor is 
included in the differential cross section formula, which depends on a Gaussian distribution of 
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the transverse momentum of the charmonium |7] . We have used Weizsacker- Williams equivalent 
photon approximation for the photon distribution of the electron [8, 9 , to calculate the cross 
section for the process e + p' — >• e + J/ip + X at low virtuality of the photon. The underlying 
partonic process at LO is jg — > cc and therefore, the only k± dependent pdf appearing is the 
gluon Sivers function. For a complete calculation of photoproduction of J/ip one has to consider 
higher order contributions and also the resolved photon contributions [7]. Also the gauge links 
or Wilson lines present in the TMD distributions are important at higher order [ID] . 

According to CEM, the cross section for charmonium production is proportional to the rate 
of production of cc pair integrated over the mass range 2m c to 2mo 

1 f2m D i 

where m c is the charm quark mass and 2m d is the DD threshold, M 2 is the squared invariant 
mass of the cc pair. 

To calculate SSA in scattering of electrons off a polarized proton target, we assume a gen- 
eralization of CEM expression by taking into account the transverse momentum dependence 
of the Weizsacker- Williams (WW) function and gluon distribution function. The numerator of 
the SSA can be written as 



d 4 ^ d 4 a± 



dydM 2 d 2 q T dydM 2 d 2 q T 



- J [d 2 k ±1 d 2 k ±g ]A N f g/pt (x g ,k±g)f 7 / e (x 1 ,k±i) 

x 5 2 {k^ + k ±g - q T )a^ cE (M 2 ) (2) 



where y is the rapidity and g-r in the transverse momentum of the charmonium; A N f g / p t (x g , k± g ) 
is the gluon Sivers function, / 7 / e (x 7 , fcj_ 7 ) is the photon distribution of the electron, given in 
the WW approximation. The denominator would have a similar expression involving the unpo- 
larized gluon distribution of the proton; f g / p (x g , k± g ), for which we use a gaussian form of k± 
distribution and a similar gaussian form for the transverse momentum dependence of the WW 
function. To extract the asymmetry produced by the Sivers function, we calculate the weighted 
asymmetry [llj 

J d<f> qT [den — da^) sm(<p qT — 4>s) 
J d<f>q T [da^ + da^} 



where <j) qT and 4>s are the azimuthal angles of the J ftp and proton spin respectively. For the 
gluon Sivers function we have used a model in our analysis, which has been used in the literature 
to calculate SSA in semi-inclusive deep inelastic scattering (SIDIS) [12] and DY process [T3] 
(see [2] for details). The parameters are taken from [T5]. Other parameters we use are 
(k 2 ±g )~(ki 1 )=Q.25GeV 2 . 

Also it is to be noted that in the model we consider for charmonium production, namely the 
color evaporation model, the only relevant scale is M 2 which is the invariant mass of the heavy 
quark pair. This is integrated between a narrow region, from Am 2 to Amjj irrespective of the 
center-of-mass energy of the experiment. So the scale evolution of the TMDs is not expected 
to affect the asymmetry too much. 

In Fig. 1 we have shown a comparison of the y and qx dependence of the asymmetry at 
JLab, HERMES, COMPASS and eRHIC. Model I refers to the parametrization in [13] and (a) 
refers to the parametrization of the gluon Sivers function in terms of an average of the u and 
d quark Sivers function jl4) . Different experiments cover different kinematical regions, and our 
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Figure 1: (Color online) The single spin asymmetry A s ^^ qT 



for the e+p^ — > e + J/ip + X 



as a function of y (left panel) and qr (right panel). The plots are for model I (a) (see text) 
compared for JLab (y/a = 4.7 GeV) (solid red line), HERMES {y/s = 7.2 GeV) (dashed green 
line), COMPASS («js = 17.33 GeV) (dotted blue line), eRHIC-1 (y/s = 31.6 GeV) (long dashed 
pink line) and eRHIC-2 = 158.1 GeV) (dot-dashed black line). 



results clearly show that the asymmetry is sizable, and that it is worthwhile to look at SSA's 
in charmonium production in order to extract information on the gluon Sivers function. 
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